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We report a first principles investigation of photocurrent generation by graphene PN junctions. 
The junctions are formed by either chemically doping with nitrogen and boron atoms, or by con- 
trolling gate voltages. Non-equilibrium Green's function (NEGF) formalism combined with density 
functional theory (DFT) is applied to calculate the photo-response function. The graphene PN 
junctions show a broad band photo-response including the terahertz range. The dependence of the 
response on the angle between the light polarization vector and the PN interface is determined. Its 
variation against photon energy Eph is calculated in the visible range. The essential properties of 
chemically doped and gate-controlled PN junctions are similar, but the former shows fingerprints of 
dopant distribution. 

PACS numbers: 



I. INTRODUCTION 

The extraordinary electronic and optical properties of 
graphene make it a promising material for novel appli- 
cations in photonics and optoelectronics.^ Experimental 
investigations have discovered several very interesting be- 
haviors of graphene - some are quite unique, such as 
the universal optical conductance, high transparencjEi, 
and broadband non-linear photoluminescenc^. Possible 
applications include transparent electrodes, photovoltaic 
cells, and touch-screen^^. One of the most interesting 
applications is the graphene photo-detector. The gap- 
less nature of graphene implies that there may always be 
an electron-hole pair in resonance with any excitation, 
allowing a broad operational wavelength ranging from 
SOOriTO to 6/im which can overcome t he long-wavelength 
limit of conventional semiconductor's^. Once excited, 
due to the high carrier mobility of graphene, hot elec- 
trons can transfer energy to the entire system in very 
short time scales (~ Ips)^, and a IQGbits^^ optical 
data link for err or-free detection of optical data stream 
has been reportecP^'^. These results demonstrate a great 
potential of graphene to be an advanced material for light 
sensors and high speed communications. 

A dipole potential is necessary for the photovoltaic 
(PV) process that produces a net dc electric current. For 
a graphene system, it can be achieved by either tuning 
a gate voltage locally or chemically doping acceptor and 
donor atoms at different regions. Photocurrent produced 
by both schemes was reported previousljli^lli^. The 
chemically doped graphene PN junction is more com- 
mon in experimental studies; the gate controlled junc- 
tion is more difficult to fabricate since it requires more 
precise setup of the device structure and light illumina- 
tion. For the chemical junctions, the doping level lies 
within hundreds of millivolts which corresponds to tera- 
hertz or far-infrared frequencies; on the other hand, the 
gate controlled junction is able to reach the higher fre- 
quency regime. 

The photo-response behavior of graphene PN junctions 



is very important and warran t syst ematical theoretical 
investigations. Recently, Refs.'^^'i^ reported theoretical 
analysis of gate controlled graphene PN junctions by solv- 
ing a model Hamiltonian and determining photo-response 
within a semi-classical transport model. To the best of 
our knowledge, photo-response of graphene PN junctions 
have not been investigated from atomistic first principles 
point of view. It is the purpose of this work to fill this 
gap and we investigate the photo-response function for 
both the chemical and gate controlled systems. In par- 
ticular, photo-response in a broad frequency range in- 
cluding the terahertz is investigated. Such a broad band 
photo-response is difficult to reach by conventional semi- 
conductors. The dependence of photo-response on the 
angle between the direction of light polarization and the 
PN interface is determined. This dependence is found to 
be well consistent with the previous theoretical resulPl. 
We also determined the dependence of photo-response 
on photon energy Ep}^ at the range of Ep^ ~ U , where 
U is the potential decline in the depletion region of the 
PN junction. In the solar visible range, for Eph < U the 
photo-response is almost linear, agreeing with that found 
in Ref.^^, and it tends to saturate when Ep^ exceeds U . 
Signatures of dopant distribution can be found in the 
photo-response function of the chemical junctions. By 
investigating both chemical doping and gate control, the 
PV properties of graphene PN junctions are compared 
and a comprehensive understanding is achieved. Finally, 
our atomistic approach provides a benchmark result for 
first-principle optoelectronic studies of the graphene sys- 
tem. 



The rest of the paper is organized as the following. In 
the next section, the calculation method is presented. 
Sections III and IV present results for the chemically 
doped and gate controlled junctions, respectively. The 
last section is reserved a for short discussion and sum- 
mary. 
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II. METHOD 

A two-probe graphene PV device consists of a scat- 
tering region sandwiched by the left and right graphene 
electrodes. Figjljplots such a two-probe PV system show- 
ing a chemically doped graphene PN junction. Here, a 
PN junction is formed by doping nitrogen (N, blue) on 
the left of the junction and boron (B, pink) on the right. 
The scattering region of the device consists of the PN 
junction and many layers of (doped) graphene on the 
left and right. Far from the PN junction, the scatter- 
ing region connects to the left and right electrodes: the 
electrode atoms are shown in the shadowed boxes. The 
electrodes extend to electron reservoirs at z = ±00 where 
photo-current are collected. In the analysis, the scatter- 
ing region should be large enough (12 atomic layers in 
our calculations) so that the electronic structure of the 
electrodes are not affected by the charge transfer at the 
PN junction in the middle of the structure. 




FIG. 1: (color online) Schematic atomic structure of the PN 
junction. The system is a graphene sheet doped with N (blue) 
atoms on the left side, and B (pink) atoms on the right. The 
regions in shadow are part of the electrodes which extend pe- 
riodically to z ± 00. DC current flows from the left to the 
right (along the z-axis) or vice versa. In the first principles 
calculations, the device is treated as periodic in the x- and 
y-directions. In the y-direction, a large vacuum region is in- 
cluded in the calculation super-cell that efi^ectively isolates 
interactions between the graphene sheet and its periodic im- 
ages. P and Q represent positions of dopant atoms at the 
interface. 

Because current flow is intrinsically a nonequilibrium 
transport problem, our calculation is based on carrying 
out real space density functional theory (DFT) calcu- 
lation within the Keldysh nonequilibrium Green's func- 
tions (NEGF) formalism^''. Very briefly, in the NEGF- 
DFT formalism, DFT is used to self-consistently calcu- 
late the electronic structure and Hamiltonian of the de- 
vice; NEGF is used to determined the quantum statisti- 
cal information that is needed to populate the electronic 
structure and calculate the density matrix; real space nu- 
merical methods are used to handle the transport and 
electro-static boundary conditions at the interface be- 
tween the electrodes and the device scattering region. 



Since the NEGF-DFT method has been well documented, 
we refer interested readers to the original literature for 
more detail^i^, in this work, a GPAW implementatioiJ^ 
of the NEGF-DFT is used for the calculations. 

Our interest in this work is to analyze the photo- 
response of the PV device. A linear polarized light is 
shined on the scattering region and not the electrodes. 
This is reasonable because very long graphene nanorib- 
bon (> lOOOnm) can be fabricated experimentallj^^ and 
the electrodes can be covered by optically nontranspar- 
ent materials. To obtain the corresponding photocurrent, 
we first calculate the Hamiltonian Hq of the two-probe 
PV device (e.g. FigjT| with out l ight, using the NEGF- 
DFT self-consistent methocP^Esl Here, the exchange- 
correlation is treated with the PBE functional^ and a 
(4,1) k-mesh is applied the in x-y plane for k-sampling. 
A single zeta polarized (SZP) atomic orbital basis set is 
used to expand all the physical quantities. For graphene, 
SZP basis is sufficient for obtaining accurate results. 

After Hq is self-consistently calculated, electron- 
photon interaction is added to it as a perturbation for 
the subsequent analysis of photocurrent. This way, the 
total Hamiltonian of the electron-photon system is: 

H = Ho + —A • p (1) 

where A is the polarization vector of the light. Assum- 
ing the electromagnetic field of a photon is a single-mode 
monochromatic plane wave, also assuming the suscep- 
tibility and dielectric constants to be homogenous, we 
haveP 

A{t) = a(^^/.)i/^(5e— * + b^e^^') (2) 

where /z^, Sr, e are the relative susceptibility, relative 
dielectric constant and absolute dielectric constant, re- 
spectively, a represents polarization of the field, /^^ is 
the photon fiux defined as the number of photons per 
unit time per unit area. 



Vy/lIrCr 

The first order electron-photon Keldysh self-energies 
arfpMSl 

^frniE) =Y.Mij,M,^[NG^^{E±hu)+{N+l)G^^{EThLo)] 

pq 

(4) 

where M,,„ = i^^i^i^^I^f'^ < l\Pz\m >, assuming 
z to be the light polarization direction. The Keldysh 
Green's functions are written as: 

G< = GO'-(*ri/i + tTnfn + S<JGoa , (5) 

G> = G"''(iri(l - f)L + tTnil - f)R + S> )Goa (6) 

where F represents the coupling of the device scattering 
region to the electrodes, namely the linewidth function 
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of the electrodes. and G"" are the retarded and 

advanced Green's functions (without photons), respec- 
tively. In Eq.(4]), G"^ is calculated by Eq.(5[6) without 
the electron-photon interaction. 

With the ground state Haniiltonian Hq of the two- 
probe PV device, quantities T, G°'', G^", G°^ can be 
calculated^^. Afterward, is calculated via Eq.(j4]), 
which can then be put into Eq.([5]) and Eq.Q to deter- 
mine the Keldysh Green's function G^. 

Finally, following ReflSOl the linear phase-coherent 
photocurrent is calculated by the following formula, 



Tr{[r^ 



r"ie)]G<{ey 



ie f de 

[/L(6)r^(6) - Jn{e)T^ie)W{e) - G'^{e)]] (7) 



Since the photon field is taken as a perturbation to Hq, 
the photocurrent obtained this way is in the small inten- 
sity regime and nonlinear optical effects are not consid- 
ered in this work. 



III. CHEMICALLY DOPED JUNCTIONS 

For the chemically doped graphene PN junction, the 
schematic atomic structure is shown in Fig|l] where a 
hetero-j unction is formed by doping the two sides with 
nitrogen or boron atoms, respectively. We investigated 
six junctions with different dopant configurations (inset 
in Fig. 6). For these junctions, while it is more realistic to 
use metal leads far away from the PN junction to collect 
current and connect the PN junction to the outside world, 
to reduce the computational demand we extended the 
doped graphene to z = ±oo to act as device leads. This 
is acceptable because photocurrent is generated in the 
PN junction where the dipole potential locates. 

To be specific, let's consider a particular configura- 
tion (the sixth pattern in the inset of Fig. 6). Fig|2] plots 
the calculated ground state properties without photons. 
The partial density of states (PDOS) in the upper panel 
of Fig|2] shows two characteristic peaks which are con- 
tributed by the donor (N) and acceptor (B) atoms just 
above and below the Fermi level (dashed vertical line). 
As a donor, the N atoms give electrons to the graphene 
which is indicated by its PDOS peak being above the 
Fermi level (blue peak near 0.6eV). Similarly, the B atoms 
accept electrons and its PDOS peak is located below the 
Fermi level (red peak near -0.9eV). On the other hand, 
a Bader charge analysis shows that the N atom actually 
gains some charge, 1.1 9e on average; and the B atom 
loses about 1.97e (see table in the upper panel of Figj2]) - 
these are in agreement with previous literature^""^. Hence 
the Bader charge analysis appears to contradict to N and 
B atoms being donors and acceptors. A further detailed 
analysis showed that the carbon atoms in the N side of 
the junction lose electrons in total but gain electrons in 
their py orbital which forms conjugate tt states for charge 
carriers that contribute to transport. On the B side of 



the junction, the carbon atoms lose tt electrons. In other 
words, the N atoms lose some high energy level charges as 
donor while attract more low energy level covalent elec- 
trons due to greater electron affinity. The situation of B 
atoms is analogous to that of N. This way, by separat- 
ing the transport carriers (the tt electrons) from the total 
electron number on the atoms, the Bader charge analy- 
sis is reconciled with the roles of N and B atoms being 
donors and acceptors. 

The carbon atoms in the system, to some extent, fol- 
low the behavior of the dopants, i.e. the C atoms in the 
N or B side form an extra conduction or valence band, re- 
spectively. Then, an electron in the extra valence bands 
(contributed by (B) and [C(B)]) can absorb a photon 
and hop to the extra conduction bands (contributed by 
(N) and [C(N)]) across the Fermi level, eventually giv- 
ing rise to the electron fiow from the B to the N side 
of the PN junction. This is also refiected in the aver- 
age effective potential shown in the lower panel of Fig|2] 
where we observe a potential decline about ^ lAeV at 
the interface (the oscillation is due to the ionic poten- 
tial). An electron-hole pair is generated after absorbing 
a photon, the electron fiows along the potential decline 
and the hole in the opposite direction. A photocurrent is 
therefore generated in the graphene PN junction by the 
photovoltaic effect without external bias. 
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FIG. 2; (color online) Upper panel: the PDOS of N and B 
atoms. There is a donor peak of N around 0.6eV (blue peak) 
and an acceptor peak of B around -0.9eV (red peak). The 
dashed vertical line at energy zero is the Fermi level. The 
result of Bader charge transfer analysis is listed in the inset 
table. Lower panel: averaged effective potential profile along 
the transport direction. Blue dashed line represents the level 
in the left electrode and red dash line for the right. 

We now turn on the light and investigate the depen- 
dence of photo-response to light polarization. Figj3]plots 
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the photo-response function. 



f = 

■' ~ eF 



(8) 



versus photon energy with different polarization direc- 
tions of the Hght. Here F is the photon flux: in our 
perturbation theory, photocurrent scales with F. 9 in 
Fig|3]is the angle between the direction of light polariza- 
tion and the y-axis. The photocurrent almost vanishes 
completely when 9 = 0, i.e., when the electric field of 
the light vibrates along the y axis. This is because the 
angular momentum quantum number of such a photon 
is equal to one along the y-axis, while both the valence 
band in the B side and the conduction band in the N 
side of the PN junction are in the tt orbital, also having 
unity angular momentum quantum number along the y- 
axis. The photon induced hoping from the valence band 
to the conduction band is therefore prohibited due to the 
mismatch of the three angular quantum numbers. This 
limitation starts to be relieved when 9 increases from zero 
and is completely removed when 9 reaches 7r/2 at which 
the maximum photocurrent is reached. As a direct result 
of the Fermi's golden rule, the photocurrent scales with 
sin 9 (inset of Figjs]) which is in good agreement with 
results of Ref, 15 Since the angle dependence is deter- 
mined by the angular momentum of the tt orbital, the 
gate controlled PN junctions to be discussed below show 
the same trend. 
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FIG. 3: (color online) Photo-response function with differ- 
ent light polarization directions. The polarization direction 
changes continuously from the tt orbital direction (y-axis) to 
the transport direction (z-axis). Inset: photo-current at pho- 
ton energy 2.0eV, versus different polarization angles. The red 
dashed line is fitting to the trigonometric function I^/2sin^9. 

Different from the discrete features predicted for car- 
bon nanotube^22]^ photo-response of graphene is smooth 
over a broad frequency range corresponding to the con- 
tinuous spectrum of the graphene material (e.g. Dirac 
bands near the Fermi level). In other words, there is 



no obvious transition features between discrete states in 
the / — Eph curve. Photons with energy Eph can excite 
electrons in the range \—Eph,Q] to generate electron-hole 
pairs; these pairs can be separated in the depletion region 
of the PN junction and contribute to the photocurrent. 

We have also investigated the role of dopant distribu- 
tion by comparing results of six samples having different 
dopant locations (only the positions indicated by letters 
P and Q in Figjljare randomly changed). These results 
will be compared with those of the gate controlled PN 
junctions in the next section. 



IV. GATE CONTROLLED JUNCTION 
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FIG. 4: (color online) (a) Gate controlled PN junction: it is 
a perfect graphene sheet with two local metal gates (shown 
as green shadowed boxes) on top. Coordinates are defined 
the same way as that in Fig[l] (b) Schematic plot of a dipole 
potential along the transport (z-axis) direction. AX z = ±oo, 
the potential equals to that of perfect graphene. Dashed lines 
indicate the slowly decaying potential over a long distance, (c) 
The DOS of gate controlled graphene PN junction. Above the 
horizontal black line is for a positive gate voltage (right side) 
where the linear Dirac bands are shifted upward; below is for 
a negative gate voltage (left side). The green vertical dashed 
line indicates the Fermi level. Shadowed areas represent the 
region where electrons can be excited by photons to contribute 
to the photocurrent. 



The structure of a gate controlled device is shown in 
Fig|4^. We imagine that a metal gate of finite width 
D is attached on either side of the graphene forming a 
transport junction. When opposite gate voltages are ap- 
plied on the two gates, a potential drop is established 
across the junction. Far left to the left-gate and far right 
to the right-gate, the system is not affected by the gate 
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potentials so that the two electrodes far away still have 
the same chemical potential, as schematically shown in 
the potential profile in Fig|4]3. Without photons, there 
is no dc current because the electrochemical potentials 
of the left and right reservoirs are the same. With pho- 
tons, electrons may be excited from the valence DOS (as 
indicated by the shadowed area in Fig|4|;) to the conduc- 
tion DOS, and a net dc current is generated by the local 
electric field. 

In practical calculations, gate voltages can be applied 
as the electrostatic boundary conditions for the Hartree 
potential^Sl For graphene, the potential decays very 
slowly due to poor screening of low dimensionality as 
schematically indicated in Fig|4j3: very long graphene 
sheet between the two gates would be needed in the cal- 
culation to correctly capture the potential profile. This 
computation difficulty can be bypassed approximately. 
In our analysis, we neglect electron-hole pair splitting in 
the decaying region of the potential, this is well justified 
due to the weak local electric field. We also neglect any 
change (due to the gate) in the self-energy: this is similar 
to the wide-band approximation. Hence, the gate volt- 
ages can be simulated in the NEGF-DFT method very 
simply as follows. We first carry out a NEGF-DFT self- 
consistent analysis of a homogeneous graphene, increas- 
ing the bias voltage to U, and then we reset the electro- 
chemical potentials of both electrodes to the Fermi level 
of the unbiased graphene in subsequent photocurrent cal- 
culations. As a result, the junction produces photocur- 
rent when the scattering region (between the gates) is 
under light illumination. 
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FIG. 5: (color online) Photo-response function of the gate 
controlled graphene PN junction under different gate voltage 
differences U displayed in the legend. On the right of this 
figure (higher photon energy range), the curves are ordered 
from low U to high [/, namely the lowest curve is 17 = 0.2V 
and highest curve is (7 = 1.4V. 

In Fig|5]we plot the photo-response function / (Eqjsj) 
subjected to different gate voltage differences U for the 



gate controlled junction, versus the photon energy E^h = 
hio. When < U, f increases linearly with Eph- In 
Refs llSlfBl by solving a semiclassical transport model the 
photo-response function / = under certain light in- 
tensity S, is also found to be linear in photon energy, 
/ = "^4^ oc u when hu ^ U {W is the width of the 
strip and f3 the slope of the potential). Here we give a 
simple argument of the linear relationship when Ep^ ~ U 
- based on the intrinsic linear dispersion of graphene. 
As shown in Fig|4|:, the "V" shaped DOS of a pristine 
graphene is shifted above or below the Fermi level by 
the gate voltage, and the distance between the two dips 
(labeled C and V) is the difference U between the gate 
voltages. A photon with energy Eph can excite valence 
electrons having energies in the range [Eph, 0] to the con- 
duction band, i.e., the number of excited electrons N^^ 
is proportional to the shadowed area indicated in Fig|4j:. 
The excited electrons are collected by electrodes at the 
edge of the scattering region where the coupling Hamil- 
tonian can be simplified as iT (F is a constant represent- 
ing the lifetime of electron at the edge). Photocurrent 
is therefore determined by two factors, Iph oc N^x^- In 
Fig|4j;, the red shadowed area means that the electron is 
excited from the right side to the left side of the junction, 
and the opposite is true for the blue shadowed area. The 
photocurrent contributed by these two areas have oppo- 
site sign, and the net current depends on their difference. 
Due to the linear dispersion of graphene, the difference 
of these two areas scales with Eph, therefore the photo- 
response - photocurrent induced by one photon, is linear 
to Eph- For Eph > U, the response involves more com- 
plicate factors rather than the n bands. From the results 
we conclude that it tends to saturate and appears quite 
flat over a broad range of Epi^. 

A very interesting issue is how does gate controlled de- 
vice compare with the chemically doped ones. To this end 
we calculated six samples of chemically doped devices in 
which the dopant locations are different (positions indi- 
cated by letters P and Q in FigjlJ. The results of the gate 
controlled device and the six chemically doped devices are 
shown in FigjHj At low energy, the photo-responses of 
all the devices are similar, suggesting that the response 
is mostly determined by the overall potential slope of 
the PN junction, namely the carbon system provides the 
main contribution to photocurrent. When photon energy 
is larger than leV, the photo-response starts to differ 
from each other. The chemically doped junctions have 
some specific undulations (black curves in Fig[6| which 
can be traced to dopant scattering. This is a character- 
istic fingerprint owing to quantum interference between 
different electron propagating paths. The photo-response 
of the gate controlled junction is smooth since the po- 
tential profile is relatively smooth"^^. It also has a larger 
response in the high frequency regime. This is because in 
chemically doped junctions, the strong local fields near 
the dopants may not contribute to photocurrent if the 
donor and acceptor sit far apart from each other and, on 
average, for chemical junctions the active carbon system 
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suffers from a weaker global field across the PN junction 
as compared to that of the gate controlled junctions. 




Photon Energy (eV) 

FIG. 6: (color online) Black solid curves are the photo- 
response function for junctions with different dopant config- 
uration (interface structures shown in the inset). These con- 
figurations are generated by changing the locations of P, Q 
sites in FigjT] The red dotted curve is the photo-response of 
the gate controlled PN junction having an average potential 
decline of 1.4eV. The polarization angle of the incident light 
is set to ^ = 7r/2 for all cases. 



V. DISCUSSION AND SUMMARY 

The calculated values of the photo-response function 
imply that for a ~ 1mm wide graphene PN junction(the 
light is shined perpendicularly to the junction plane), 
at the standard solar light intensity of Q.lW/cm? and 
a photon energy of leV , the corresponding photocurrent 
should be in the order of nA which is detectable exper- 
imentally. For the chemically doped graphene PN junc- 
tions with high dopant concentration, short-range disor- 
der may trigger charge localization, coherent back scat- 



tering and damped quantum interference^. The localiza- 
tion length at the Dirac point of graphene is estimated to 
be ^ 200nm by two-dimensional scaling theorji^D. The 
predicted photocurrent may be suppressed to some ex- 
tent by such effects. 

We mention in passing that effects of charge-transfer 
exciton which is important for organic PV devices has not 
been considered in our analysis because it is negligible for 
single layer graphene, especially in the low photon energy 
regime (< leV)^-. Finally, we did not observe quantum 
interference between two electron paths accompanied by 
resonant absorption of photons discussed theoretically in 
Ref(271 The reason may lie in the mismatch of the de- 
vice parameters which are required to satisfy both strong 
reflection and smooth potential profile, hence to observe 
that phenomenon it requires a device size in the order of 
lOOnm which is much larger than the systems we inves- 
tigated here. 

To summarize, we have carried out a first principles 
atomistic analysis of the photo-response of graphene PN 
junctions. For both chemically doped and gate con- 
trolled junctions, we examined several main features of 
the photo-response function. Our calculations show that 
the graphene PN junctions have a broad band photo- 
response - including terahertz, which is difficult to reach 
by conventional semiconductors. The dependence of 
photo-response on the angle between the direction of 
light polarization and the PN interface is consistent with 
the previous theoretical predictiorP^. The overall trend 
of the dependence of photo-response against photon en- 
ergy, except the chemical fingerprints observed in the 
chemically doped junctions, agrees well with that of the 
semiclassical model in the small intensity regime. Impor- 
tantly, the essential properties of the photo-response for 
the two kinds of PN junctions are found to be similar. 
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